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It  is  shown  in  this  study  that  the  relationships  of  the  delay 
times  of  secondary  arrivals  to  the  distance  and  azimuth  of 
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ferentiating the  source  characteristics  of  an  event.  Cepstrum 
analysis  is  the  method  used  to  detect  the  delay  times  of 
secondary  arrivals  and  was  applied  in  a first  zone  study  to  I 
■®l«cted  Nevada  Test  Site  events  and  one  natural  earthquake.  1 
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\ 20.  Abstract  (cont'd) 

"^The  theoretical  basis  of  cepstrum  analysis  was  extended  to 

include  three  secondary  arrivals,  and  a formula  was  establishefd 
for  recovering  the  scaling  factor  of  secondary  arrivals  with 
respect  to  the  initial  P arrival  from  an  event. 

The  results  indicate  that  cepstrum  analysis  can  be  used  to 
detect  the  delay  time  of  secondary  arrivals,  but  it  is 
difficult  to  interpret  the  delay-distance  and  delay -azimuth 
relationships  and  thereby  identify  the  secondary  arrivals. 

The  detection  of  multiple  events  appeared  possible  in  some 
cases,  but  station  coverage  requirements  are  severe. 

Knowledge  of  within  5 00  km  of  a source  has  possibly 
been  extended.  ^ 
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Technical  Report  Summary 

This  is  the  third  semiannual  report  dealing  with  an 
investigation  of  multiple  seismic  events  and  first  zone 
discriminants.  Reported  here  are  the  results  obtained 
from  studying  secondary  arrivals  from  a series  of  under- 
ground nuclear  shots  at  the  Nevada  Test  Site  and  a natural 
earthquake  in  the  same  general  source  region  utilizing 
Cepstrum  analysis  techniques. 

Cepstrum  analysis,  when  applied  to  real  data  at  first 
zone  distances,  yields  a number  of  peaks  in  the  resulting 
cepstra.  The  signal  to  noise  level  appears  rather  low  and 
the  peaks  cannot  be  readily  differentiated  on  the  basis 
of  amplitude.  Lacking  any  significant  differences  in 
amplitude,  it  is  necessary  to  consider  all  of  the  resulting 
peaks.  When  the  delays  associated  with  all  the  peaks  are 
considered  as  a function  of  station  distance  and  azimuth, 
tentative  identifications  are  possible.  The  scatter  of 
the  resulting  data,  however,  is  high  and  some  subjective 
evaluation  is  required.  In  view  of  the  subjective  ele- 
ments in  the  identification,  it  does  not  appear  that  the 
method  can  be  converted  to  a simple  machine  procedure. 

The  identification  of  arrivals  does  appear  possible 
vdien  sufficient  data  are  available  and  subjective  evalua- 
tions are  allowed.  Selected  arrivals  exhibit  the 
characteristics  of  single  sources  and  multiples  of  hori- 
zontally separated  events. 
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A delay  arrival,  occurring  at  approximately  0.70  sec 
at  200  km  with  moveout  of  about  0.08  sec  per  100  km, 
occurred  for  five  of  the  six  nuclear  events.  This  is 
tentatively  identified  as  pP.  Very  little  data  is  avail- 
I able  on  pP  at  these  distance  and  depth  ranges.  Thus,  if 

this  identification  is  correct,  the  knowledge  of  pP  is 
extended. 

Due  to  the  necessity  of  considering  the  measured 
delays  as  a function  of  distance  and  azimuth,  station 
coverage  requirements  for  the  detection  of  multiple  events 
are  severe.  The  detection  of  P-P  delays  appears  to  require 
stations  at  20®  azimuth  intervals.  The  detection  of  the 
pP  arrival  requires  at  least  8 stations  located  at  50-75 
km  intervals  over  a 500  km  range. 

The  "a"  scaling  factor  was  not  recovered  for  the 
nuclear  events  considered  in  this  study.  The  failure  to 
recover  the  value  resulted  from  limitations  in  available 
computer  time.  In  addition,  the  available  station  coverage 
of  azimuth  and  distance  was  somewhat  limited  and  far  from 
the  ideal  given  above.  Despite  these  limitations,  the 
^ detection  of  multiple  events  appeared  possible  in  some 

cases . 

On  the  basis  of  theoretical  considerations  and  limited 
tests,  values  of  the  ”a”  factor  can  be  recovered  with  some 
degree  of  accuracy.  While  recovery  is  possible,  the 
resulting  "a"  is  both  theoretically  and  practically 


ambiguous,  and  because  of  this  may  prove  of  little  value 
in  the  identification  of  multiple  events.  However,  despite 
the  ambiguous  nature  of  "a",  this  value  would  be  useful  in 
the  identification  of  pP. 

No  obvious  propagation  path  effects  due  to  differences 
in  geology  are  evident  from  the  data.  While  the  results 
from  an  actual  earthquake  appear  different  from  the  results 
obtained  from  nuclear  events,  insufficient  data  was  avail- 
able to  establish  a quantitative  difference.  It  is 
interesting  to  note,  however,  that  for  the  available  data, 
the  variance  associated  with  the  predominant  peaks  of  the 
earthquake  was  0.5  seconds,  while  the  variance  associated 
with  the  predominant  peaks  of  nuclear  events  was  approxi- 
mately 0.1  seconds. 
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Dear  Sir: 

This  report  is  a summary  of  research  dealing  with  multiple 
seismic  events  and  first  zone  discriminants.  The  research 
is  divided  into  the  following  categories  and  will  be  dis- 
cussed individually. 
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Introduction 


The  purpose  of  the  research  reported  here  is  to 
evaluate  one  of  the  methods  for  the  discrimination  of 
multiple  nuclear  events  from  single  nuclear  and  natural 
events  at  distances  less  than  1000  km.  The  basis  for  the 
method  of  discrimination  investigated  is  the  relationship 
of  the  delay  of  arrivals,  following  the  initial  P arrival, 
to  azimuth  and  distance.  These  relationships  should  be 
characteristic  of  the  event  type.  Cepstrum  analysis  is 
used  to  determine  the  delay  of  the  secondary  arrivals 
(Bogart,  Healy,  and  Tukey,  1973).  This  method  of  analysis 
will  be  applied  to  records  of  selected  events  within  the 
Nevada  Test  Site.  This  study  will  be  a near  source  appli- 
cation of  cepstrum  analysis  to  events  at  distances  ranging 
from  150  km  to  415  km.  Where  possible,  the  relationships 
of  any  detected  arrivals  to  distance  and  azimuth  will  be 
used  to  predict  the  nature  of  the  event. 

This  report  is  based  on  the  thesis  prepared  by 
Zimdars  (1974)  which  was  sponsored  by  this  project.  The 
results  of  our  investigations  into  spectral  and  correlation 
analyses  between  the  seismic  signals  from  earthquakes  and 
nuclear  events  will  be  presented  in  the  next  scheduled 
semiannual  report. 


Theoretical  Considerations 

Source  Parameters  and  Multiple  Arrivals 
One  method  of  determining  the  multiple  nature  of  an 
event  is  to  detect  the  multiple  nature  of  the  direct  P 
arrivals  which  occur  from  the  multiple  sources.  Associated 
with  a multiple  source,  there  is,  thus,  a P-P  delay  time, 
the  time  difference  between  the  arrival  of  the  direct  P 
wave  from  one  event  and  the  arrival  of  the  direct  P from 
a second  event-  The  major  problems  of  this  approach  are 
caused  by  the  existence  of  contaminating  secondary  arrivals. 
For  natural  earthquakes  and  single  nuclear  sources  there 
exists  a surface  reflected  P wave,  pP,  of  similar  shape  but 
different  polarity.  In  addition  to  pP  for  a single  nuclear 
source,  there  is  another  possible  secondary  arrival  caused 
by  the  slapdown  of  surface  materials  (King  and  others, 

1974;  Springer,  1974).  Thus,  each  single  nuclear  explosive 
source  may  consist  of  three  P phases,  a direct  P,  pP,  and 
a slapdown  arrival  Pg,  and  these  have  the  associated  delay 
time,  pP-P  and  Pg-P.  All  of  these  arrivals  will  occur 

within  the  normal  P wave  coda. 

A multiple  source  composed  of  two  events  consists  of 
two  sets  of  the  three  P phases  and  has  the  associated  delay 

times,  pP;^^-Pj^,  P^2"^l'  ^S2~^l* 

magnitude  of  the  delay  time  and  the  order  of  the  arrival 

of  the  phases  at  a particular  station  is  dependent  upon 
the  location  of  the  events  and  the  shot  times. 
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The  possibility  of  the  existence  of  a number  of 

\ 

secondary  arrivals  makes  it  difficult  to  relate  a secondary 
arrival  from  a single  record  to  any  source  characteristics. 
The  different  secondary  arrivals  must  be  differentiated 
from  each  other.  A possible  method  of  differentiating 
the  secondary  arrivals  exists  in  their  differing  delay  to 
distance  and  delay  to  azimuth  relationships.  Thus,  the 
use  of  the  arrivals  from  a number  of  stations  located  at 
different  distances  and  azimuths  provide  a possible  rela- 
tionship  to  the  source  characteristics. 

In  the  following  sections  the  relation  of  delay  to 
azimuth  and  distance  will  be  considered  for  each  of  the 
various  potential  secondary  arrivals.  In  general,  to  keep 
the  involved  parameters  reasonable,  multiple  events  will 
be  considered  as  vertically  or  horizontally  separated  and 
simultaneously  detonated.  The  generalization  to  inclined 
separation  or  delayed  firing  will  normally  present  little 
•difficulty. 

P-P 

The  delay  dista.nce  curves  for  two  events  separated 
horizontally  and  two  events  separated  vertically  may  be 
quantitatively  determined  using  the  simplified  diagrams  of 
Figures  1 and  2.  Choosing  a 1 km  distance  between  events 
and  a depth  of  burial  of  .5  km,  the  calculated  distance 
delay  curves  for  three  velocities  are  shown  in  Figures  3 
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O - Midpoint  of  epicenter  line 
C - Station  on  surface  of  the  earth 
A,  B - Explosion  points 
X - Distance  to  station 
d - Depth  of  burial 
2s  - Distance  between  events 
Pi'  P2  “ of  P wave  of  velocity,  V 

= (d^  + (X  + P2  - (d^  + (X  - 


Delay  time  = 

V 


(d^  + (x  + - (d^  + (x  - 


V 


Figure  1.  Schematic  diagram  of  P wave  paths  for  two 

horizontally  separated  events.  Calculations 
are  for  the  delay  time  of  the  direct  P 
arrivals. 
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O - Epicenter 


C - Station  on  surface  of  the  earth 


A,  B - Explosion  points 


X - Distance  to  station 


d - Depth  of  burial 


s - Distance  between  events 


Pj^,  P2  - Path  of  P wave  of  velocity,  V 
p.  = (x^  + (d  + p-  ■ (x^  + 


Delay  time 


= Pi  - P2 


f r 


(x^  + (d  + - (X"  + 


Figure  2.  Schematic  diagram  of  P wave  paths  for  two 

vertically  separated  events.  Calculations  are 
for  the  delay  time  of  the  direct  P arrivals. 
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and  A.  For  two  evonts  separated  horizontally  the  delay 

approaches  constant  values  of  .40,  .25,  and  .17  seconds  i 

within  5 km  of  the  source  for  velocities  of  2.50,  4.00, 

and  6.00  km/sec,  respectively.  For  two  events  vertically 

separated  by  1 km,  the  delay  approach.es  zero  v/ithin  50  km 

of  the  source.  The  maximum  delay  in  this  case  occurs  at 

the  epicenter  of  the  events  and  for  this  model  has  values 

of  .40,  .25,  and  .17  seconds  for  velocities  of  2.50,  4.00, 

and  6.00  km,  respectively. 

In  order  to  generalize  the  delay  distance  plots, 
dimensionless  plots  were  prepared.  The  abscissa  of  the 
plot  is  the  ratio  of  the  station  distance  to  the  distance 
betv/een  sources,  and  the  ordinate  is  the  percent  of 
maximum  delay.  Dimensionless  plots  for  horizontally  and 
vertically  separated  multiple  sources  are  given  in  Figure 
5.  It  is  evident  that  the  delay  approaches  a constant 
value  within  5 source  separations  of  the  epicenter  for 
horizontally  separated  events.  For  vertically  separated 
events  the  delay  approaches  a constant  value  within  50 
source  separations  of  the  epicenter. 

For  two  events  vertically  and  horizontally  separated  j 

i 

by  a similar  amount  the  distance-delay  plot  has  the  form  | 

of  the  plot  for  horizontally  separated  events . Only  the 
beginning  section  of  the  plot  indicates  that  the  events 
have  a vertical  sepciration.  A plot  for  events  having  a 
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DELAY  (SEC) 


Figure  3.  Three  delay-distance  plots  for  horizontally 
fe"  separated  events.  Distance  between  events 

I is  1 km.  Depth  of  burial  is  .5  km.  Velocities 

"S  4%  (V)  are  2.50,  4.00,  and  6.00  km/sec. 

'’’i,  V'.‘ 
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DISTANCE  (KM) 


Figure  4.  Three  delay-distance  plots  for  vertically 

separated  events.  Distance  between  events  is 
1 kin.  Depth  of  burial  is  .5  km.  Velocities 
(V)  are  2.50,  4.00,  and  6.00  km/sec. 
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STATION  DISTANCE/DISTANCE  BETWEEN  EVENTS 
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Figure  5.  a)  Dimensionless  plot  for  two  events  separated 
horizontally. 

b)  Dimensionless  plot  for  two  events  separated 
vertically. 
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horizontal  separation  of  l o km  i ^ 

J..U  j^n)  and  a vertT/'.-ii  .. 

1.0  km  is  vertical  separation 

in  Figure  6. 

Azimuth-delay  plots  for*  -■  v, 

for  a horizontal  seoaration  of 
events  were  determined  for  di.f,  . ‘ 

5 and  1 n V -J^otances  between  events  of  1 

S-  and  10  to  and  a velocity  of  2.50  tan/seo  r d- 

and  plots  are  .K  ™/sec.  Tne  diagram 

Plots  are  shown  in  Figures  7 and  8 The  d», 
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dependent  upon  the' distance  between  sources  and  not 
the  azimuth  Tha  ^ °t  upon 

expected  delay  azimuth  plot  for  any 


50 

DISTANCE  (KM) 


Figure  6.  Delay-distance  plot  for  an  event  separated 
vertically  and  horizontally.  Vertical 
separation  is  1 km  and  horizontal  separation 
is  1 km.  Velocity  is  2.50  km/sec. 
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(-3,0, d)  (s,o,d) 


0 - Midpoint  of  epicenter  line 
A,  B - Explosion  points 

C - Station  on  surface  of  earth 
r - Distance  to  station 
d - Depth  of  burial 

2s  - Distcince  between  events 

0 - Azimuth 
V - velocity 

— 2 2 2 2 1/2 

AC  = ((rsinO  + s)  + r cos  0 + d ) 

— 2 2 9 2 1/2 

BC  = ((rsin0  - s)  + r cos^0  + d ) 

|ac  - ic| 

Delay  time  ■ — 

V 

|(r^  + s^  + d^  + 2rssin0)^^^  - (r^  + s^  + d^  - 2rssin0)^^^| 

■ ■■■“■'  " ■ ■■  - rn..™.,...  ,r..  , .r  ....... ..  l'  I I -'l 
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Figure  7.  Diagram  of  distance,  depth,  and  azimuth 

relationships  for  two  horizontally  separated 
events.  Calculations  are  for  the  delay  time 
of  the  direct  P arrivals. 
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Figure  8.  Three  delay-azimuth  plots  for  horizontally 

separated  events.  Depth  of  burial  is  .5  km. 
Velocity  is  2.50  km/sec.  Distances  between 
events  (s)  are  1,  5,  and  10  km. 
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0 - Epicenter 
A,  B - Explosion  points 

C - Station  on  surface  of  earth 
r - Distance  to  station 
d - Depth  of  burial 
s - Distance  between  events 
0 - Azimuth 
V - velocity 


AC  - (r^sin^e  + r^cos^e  + 
ic  = (r^sin^e  r^cos^f  (s  + 


Delay  time 


(r-^  + d*') 


2\V2  _ (j.2  + (3  + 


2.1/2, 


Figure  9.  Diagram  of  distance,  depth,  and  azimuth  rela- 
tionships for  two  vertically  separated  events 
Calculations  are  for  the  delay  time  of  ihe 
direct  P arrivals. 
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reasonable  vertical  separation  of  sources  is  a point  at 
the  origin . 
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pP-? 

The  pP  phase  of  a shallow  event  is  difficult  or 
impossible  to  identify  on  a seismogram.  The  1968  Seismo- 
logical  Tables  for  P phases  lists  h pP-P  interval  of  3.4 
seconds  for  an  event  at  distances  of  2“  to  13®  and  a depth 
of  15  kilometers.  An  increasing  pP-P  interval  with 
increasing  epicentral  distance  is  indicated  by  the  table 
listings.  The  station  distances  in  this  study  range  from 
1®  to  4®  and  the  shot  depths  from  .15  to  .56  kilometer. 

The  pP-P  delay  times  for  thirty-seven  U.S.  underground 
tests  have  been  estimated  by  Springer  (1974)  from  surface- 
zero  accelerograms  and  range  from  .17  to  1.00  seconds. 
Table  1 is  a listing  of  Springer's  estimates.  The  three 
lowest  estimates  of  .17/  .21,  and  .24  seconds  were  for 
shots  detonated  in  granite.  The  shot  mediums  in  this 
study  are  alluvium  and  tuff,  and  the  lowest  estimates  by 
Springer  for  shots  in  alluvium  and  tuff  are  .34  and  .33 
seconds,  respectively. 

Based  upon  the  above  information, ' the  pP-P  delay 
times  for  the  shots  in  this  study  are  expected  to  be  with- 
in the  .20  to  1.5  second  range.  The  pP-P  delay  may  be 
constant  over  the  1*  to  4®  epicentral  distance  range  or 
an  increasing  pP-P  delay  with  increasing  distance  is 


li 


i 
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possible.  IVhen  the  pP-P  delay  is  constant  over  the  epi- 
ccVitral  distance  range,  the  expected  delay-distance  plot 
' for  pP-P  is  a horizontal  line,  and  the  delay-azimuth  plot 

* is  a circle  v;ith  center  at  the  origin.  For  the  case  of 

j pP-P  increasing  with  distance  from  the  station,  the  delay- 

i 

distance  plot  is  a curve  shov/ing  the  increasing  delay  with 
distance,  and  the  delay-azimuth  plot  is  a circle  with 
center  at  the  origin  for  those  stations  of  similar  epi- 
central  distance. 


Pg-P 

Surface  material  above  an  underground  nuclear  explo- 
sion is  thrown  upward  by  spallation,  which  is  the  parting 
of  near-surface  layers  originally  in  contact  (Eisler  and 
Chilton,  1964)  . The  parting  is  described  by  Springer 
(1974)  as  occurring  when 

interference  between  the  reflected  tension  wave 
and  the  incident  compressive  stress  wave  may 
produce  a net  tensile  stress  equal  to  the  sum 
of  lithostatic  pressure  (overburden)  and  tensile 
strength  of  the  rock.  Separation  of  the  rock 
layers  above  this  point  occurs  because  some  of 
the  stress  wave  energy  is  imparted  to  tlie 
separated  layers  in  the  form  of  kinetic  energy. 

That  is,  a layer  of  the  earth  separates  and 
continues  ballistically  upward.  The  material 
eventually  falls  under  the  force  of  gravity  and 
lands  with  a considerable  shock. 

The  forcible  landing  of  the  spalled  layers  causes  the 

• slapdowu  phase.  The  occurrence  of  spallation  depends  upon 

the  yield  of  the  explosion  and  the  depth  of  burial. 


Spallation  would  not  occur  for  a vc,ry  low  yield  explosion 
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or  for  an  explosion  buried  at  great  depth.  Also,  no 
spallation  would  occur  if  an  explosion  was  not  contained. 

The  slapdown  phase  follows  the  P phase  by  the  same 
interval  at  all  stations.  The  Pg-P  delay  distance  plot 
is  a horizontal  line/  and  the  delay  azimuth  plot  is  a 
circle  with  center  at  the  origin. 

Pg-P  delays  of  .87,  1.35,  and  1.85  seconds  were 
determined  for  the  LONGSHOT,  MILROW,  and  CANNIKIN  nuclear 
explosions  by  King  and  others  (1974)  using  a spectral 
ratio  method.  Springer  (1974)  estimated  Ps~P  delays  from 
surface- zero  accelerograms  for  twenty-seven  U.S.  under- 
ground explosions.  His  slapdown  estimates  range  from  .75 
to  3.07  seconds  and  are  listed  in  Table  1. 

The  differentiation  of  the  secondary  arrivals  by 
their  delay  to  distance  and  delay  to  azimuth  relationships 
provide  a possible  method  6f  determining  the  source 
characteristics  of  an  event.  The  delay  to  distance  and 
azimuth  relationships  for  P,  pP/  and  Pg  are  svimmarized  in 
Figilre  10.  It  is  evident  from  the  figure  that  when  taken 
together  the  relationship  of  delay  to  distance  and  azimuth 
define  the  source  type  involved.  Thus,  it  may  be  possible 
to  use  plots  of  this  nature  to  differentiate  the  various 
source  types  if  the  interference  of  the  P-coda  is  signifi- 
cantly low. 

.7 
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TABLE  1 

Estimated  Delay  Times  Determined  from 
Surface-Zero  Accelerograms 
(Modified  from  Springer,  1974,  Table  1,  p.  586; 


Springer  and 
Shot  Name 

Kinnaman,  1971 
Date 

, Table  1,  p. 

Depth  of 
Burial 
(km) 

1074) 

pP-P 

(sec) 

P3-? 

(sec) 

It^INIER 

9/19/57 

.25 

0.33 

0.75 

GNOME 

12/10/61 

.36 

0.36 

1.32 

HARD  HAT 

2/15/62 

.29 

0.17 

_ ‘k 

SHOAL 

10/26/63 

.37 

0.21 

0.90 

S.^JLMON 

10/22/64 

.83 

0.58 

_ k 

11/05/64 

.40 

0.37 

0.92 

LONG  SHOT 

10/29/65 

.70 

0.45 

1.76 

DU RYE A 

4/14/66 

.54 

0.59 

1.39 

PILE  DRIVER 

6/02/66 

.46 

0.24 

1.32 

VULCAN 

6/25/66 

.32 

0.41 

HALPLE.AX 

6/30/66 

.82 

0.72 

2.93 

NEW  POINT 

12/13/66 

.24 

0.35 

GPGELEY 

12/20/66 

1.22 

0.97 

AGILE 

2/23/67 

.73 

0.83 

COM-IODORE  ■ 

5/20/67 

.75 

0.84 

SCOTCH 

5/23/67 

.98 

0.91 

1.73 

LANPHER 

10/18/67 

.71 

0.86 

G.AS  BUGGY 

12/10/67 

1.29  • 

0.87 

1.17 

HUPM03ILE 

1/18/68 

.25 

0.34 

0.99 

FAULTLESS 

1/19/68 

.98 

0.80 

2.17 

KNOX 

2/21/68 

.65 

0.67 

2.47 

30XCA.R 

4/26/68 

1.16 

0.96 

3.07 

NOGGIN 

9/06/68 

.58 

0.75 

2.92 
or  3.72 

HUTCH 

7/16/69 

.55 

0.59 

1.92 

J 0KOi.’*l 

9/16/69 

1.16 

0.98 

2.78 

MIL ROW 

10/02/69 

1.22 

0.74 

2.49 

CALABASH 

10/29/69 

.63 

0.70 

HA.NDLEY 

3/26/70 

1.21 

0.95 

2.84 

CARPETBAG 

12/17/70 

. 66 

0.67 

2.31 

HAR3ELL 

6/24/71 

.52 

0.63 

1.03 

CATHAY 

10/08/71 

.38 

0.54 

1.02 

CANN I KIN 

11/06/71 

1.79 

1.00 

_* 

CI-IAENACTIS 

12/14/71 

.33 

0.41 

0.80 

MONERO 

5/19/72 

.54 

0.69 

l.OS 
or  1.25 

DIAI40ND  SKULLS 

7/20/72 

.42 

0.47 

1.59 

OSCURO 

9/21/72 

.56 

0.69 

2.12 

DELPHINIUM 

9/26/72 

.30 

0.40 

0.90 

* Data  unavailable  or  nonexistent. 


Deiay-Azimuth 


Delay --Distance 

Delay , T 
One  event 

T = pP-P  T 


or 


T 


Two  events 
separated 
horizontally 
T = P-P 


Two  events 
separated 
vertically 
T = P-P 


Figure  10,  Summary  of  delay-distance  and  delay-azimuth 
relationships  for  P,  pP#  and  Pg. 
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Method  for  the  Detection  of  Multiple  Arrivals 

Cepstrum  analysis  was  introduced  by  Bogert,  Hcaly, 
and  Tukoy  (1963)  as  a method  for  detecting  the  existence 
and  timing  of  echoes.  The  authors  defined  the  cepstrum 
as  the  spectrum  of  the  spectrum  of  a time  series.  The 
method  involves  the  application  of  two  Fourier  transforms 
to  l;lin  daLa  nut. 

The  authors  demonstrated  the  potential  of  applying 
time  series  techniques  of  spectral  analysis  to  a series  in 
the  frequency  domain.  Using  an  artificial  time  series  with 
cin  echo,  the  echo  was  detected  by  obtaining  the  Fourier 
amplitude  spectrum  of  the  series,  determining  the  logarithm 
of  the  spectrum,  and  obtaining  the  Fourier  amplitude 
spectrum  of  the  logarithm  of  the  spectrum.  The  logarithm 
of  the  spectrum  is  a frequency  series  and  the  Fourier 
transform  of  the  frequency  series  yields  a series  in  the 
time  domain.  Thir.  is  essentially  the  spectrum  of  a spec- 
trum and  was  termed  the  cepstrum.  A peak  in  the  cepstrum 
indicates  the  existence  of  an  echo  at  the  corresponding 
time  in  the  original  series. 

For  purposes  of  completeness  and  understanding,  the 
original  work  of  Bogert,  Healy,  and  Tukey  (1963)  is  out- 
lined here.  Consider  a time  series,  f(t),  with  echo  to 
consist  of  s(t)  and  s(t)  multiplied  by  a constant  value 
a and  delayed  by  time  x.  The  series  with  echo  is  given 
by 


•L‘* 

‘A  V 
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f (t)  = s(t)  + as(t  - T)  . <1) 

ThG  Fourier  transform/  f,  of  f(t)  is 
F(m)  = S(ai)  (1  + ae^“^) 
with  the  amplitude  spectrtim  being 

2 1/2 

1F(oj)[  = lS(m)l(l  + 2acosaiT  + a ) . (3) 

Taking  the  logarithm  of  the  spectrum  converts  the  multi- 
plicative effect  of  the  echo  into  an  additive  ei.fect  and 

yields 

loglF(m)l  = loglS(u)l  + 1/2  log (1  + 2acosuJT  + a^)  . (4) 

Applying  the  log  expansion/ 

logd  + X)  = X - + — where  -1<X<1  (5) 

and  assuming  a<l  the  expression  may  be  reduced  to 

loglF(ai)l  =doglS(o))l  + acosuT.  (6) 

The  log  spectrum  is  a frequency  series.  The  Fourier  trans- 
form of  the  log  spectrum/  which  converts  the  frequency 
series  into  the  time  domain/  is  given  by 


C(t)  = FloglS(o))l  + a6(t  - t)  . 
The  amplitude  spectrxim  of  C(t)  , 


C(t)  I ■=■  I Flog  1 S (w)  I ! + a|6(t  - x) 


(8) 


V 

iG  the  cepstrura.  The  Fourier  transform  of  the  function 
cosuT  is  the  delta  function/  6 (t  - x) , which  has  an 
amplitude  of  zero  at  all  t except  t = t / v/here  its  ampli- 
tude tends  to  infinity.  Thus,  a peak  occurs  in  the 
cepstrum  at  time  t = t,  as  a result  of  the  echo  in  the 
original  series  which  was  delayed  by  time  t. 

The  condition  requiring  a to  be  less  than  1 was 
satisfied  when  applying  cepstrum  analysis  to  the  detection 
of  echoes,  since  the  amplitude  of  an  echo  is  normally  less 
than  the  amplitude  of  the  first  arrival.  The  a<l  condi- 
tion is  also  valid  for  the  detection  of  reflected  phases 
from  an  earthquake  or  single  nuclear  event  for  the  same 
reason.  An  a value  greater  than  1 is  possible  when  apply- 
ing cepstrura  ’analysis  to  the  detection  of  a multiple  nuclear 
event.  The  second  direct  P arrival  from  a multiple  event*  • 
may  have  an  amplitude  on  the  order  of,  or  greater  than, 

the  amplitude  of  the  first  direct  P. 

« 

A logarithm  expansion  which  allows  a to  be  greater 


than  1 is 
logX  = logb 


+ (X  - b)  _ (X  - b)^  ^ (X  - b)^  

b 2b2  3b3 

where  0<X<.2b. 


(9) 


2 2 

Letting  x **  1 + 2acoswT  + a and  b = 1 a and  applying 
this  expansion  when  taking  the  logarithm  of  the  spectra 
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and  disregarding  second  order  terms#  the  expression  for 
the  log  spectrvim  becomes 

loglF(to)]  = loglS(u))l  + 1/2  log(l  + (10) 

a 

+ r coswT.. 

1 + 


The  Fourier  transform  of  the  log  spectrum  is 


C(t)  = F(loglS(w)l  + 1/2  logd  + a )) 


(11) 


1 + a- 


5(t  - T) 


and  the  amplitude  spectrum#  which  is  the  cepstrum  is 


,C(t)  1 = lF(loglS(oj)  1 + 1/2  logd  + a^))  1 


(12) 


1 + a‘ 


6(t  - T)  1. 


Except  for  constant  terms  and  scaling  factors  this  is 
equivalent  to  the  original  expression  and  would  produce  an 
amplitude  peak  at  time  t.  Cepstrum  analysis  can  therefore 
be  expected  to  detect  a secondary  arrival  stronger  than 
the  first  arrival. 

While  cepstrum  analysis  has  received  some  application 
in  the  simple  detection  of  secondary  arrivals#  it  does  not 
appear  that  anyone  has  attempted  the  recovery  of  the  a 
scaling  factor  (Plinn  and  others#  1973;  Cohen,  1970). 
Knowledge  of  a would  be  useful  in  differentiating  the 
secondary  arrivals.  King  and  others  (1974)  estimated  the 
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air.plitudG  of  the  pP  and  P3  phases  of  the  LONCSIIOT,  MILRO/7, 
and  CANNIKIN  explosions  to  be  .3-. 5 times  the  corresponding 
direct  P amplitudes.  If  the  a factor  could  be  recovered 
for  the  secondary  arrivals,  any  values  of  a^.1  would 
clearly  suggest  the  arrival  of  a second  direct  P from  a 
multiple  source. 

From  the  expression  for  the  cepstrum,  equation 
the  value  of  the  peak  in  the  cepstrum  at  time  t = t 

P = lC(t)  1 - iFdog.jS  (oj)  1 + 1/2  log(l  + a^))l 


(12)  , 
is 

(13) 


1 + a' 


-l6(t  - T)l. 


In  digital  applications  the  amplitude  of  6 (t  - t)  has  a 
finite  value  which  may  be  denoted  1 6 (t  - x)  jc.  Letting 


1-+  a‘ 


(14) 


then 

Qa^-a  + Q=0, 
and  solving  for  a yields 

1 + A - 4q2 

a = • 

2Q 


(15) 


(16) 


Designating  the  roots  of  the  equation  as  a^^  and  a2,  it  can 
be  shown  that  ai  “ -i-.  The  reciprocal  property  of  the 
roots  limits  the  usefulness  of  the  recovered  a values. 
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A value  of  a in  the  range  .8<a<1.25  may  aid  in  differentia- 
tion of  a second  direct  P from  other  secondary  arrivals. 
.Values  of  a outside  of  this  range  will  be  useful  only  if 
other  information  about  the  secondary  arrivals  is  avail- 
able. 

The  theoretical  basis  for  cepstrum  analysis  and  the 
recovery  of  a discussed  above  involve  only  a single 
secondary  arrival.  To  extend  the  theoretical  base  to 
include  three  secondary  arrivals  consider  a signal 

f(t)  = s(t)  + a,s(t  - Xi)  + a2s(t  - T2)  + a3S(t  - T3)  (17) 

where  Tx'^‘^2'^'^3* 

The  Fourier  transform  is ' 

F(w)  = S (u)  (1  + .+ 

and  the  amplitude  spectrum  is 

1f(w)  I = lS(u)  1 (1  + + a2^  + 33^  + 2a2^cosuTj^ 

+ 2a2COSUT2  + 2a3COS(0T3  + 2a2_a2COsu - T2)  (19) 

1/2 

+ 2a2_a3C0su  (t^  - T^)  + 2a2a3COsu(i2  “ T3) ) . 

The  expression  for  the  log  spectrum  is 

loglF(u)l  = loglS(w)l  + 1/2  log(l  + a^^  + 33^  '*‘'^3^ 

+ 2aj_coswT2_  + 2a2COSWT2  + 2a3COSUT3  +*’2a2_a2COSW (x^  “*  13) 

+ r “*  "*^3)  2a2a3COSw(T2  - X3)). 


Applying  the  log  expanaion. 


\ 

logl  F(o))|  = ioglS(w)l  + 1/2  log(l  + a-^  + + ^3^) 

1 

+ 2 2 2 (^iCOswT]_  + fl2COsa)T2  + a^coswT^ 

1 "i"  a + 32  "t"  a^ 

+ aia2COSO)(Ti  - T2)  + 32^330030) - T3)  (21) 

+ a2a3COSw(X2  - 


yielding  as  the  cepstrum 


IC(t)I  = lF(loglS(oj)  1 + 1/2  logd  + + 32^  + 33^)! 

r+  a'2'+\  2 + ^ • ''1)1  * - '2>1 

1 + 31+32+33 

+ 33! 6 (t  - T3) I + 3^321^ (t  - (T^  - T2) ) I 
+ 33^33!  6 (t  - (T^_  - T3))  1 + 3233)6  (t  - (X2  - T3))  1 . 


(22) 


It  is  evident  that  the  cepstrum  will  become  more 
complex  with  the  addition  of  more  secondary  arrivals.  The 
peaks  occurring  at.  ^3  - x^^r  ” T^lr  '^3  ” 2 

interpretation  of  the  cepstrum  more  difficult.  This  will 
be  true  especially  in  the  case  of  a multiple  event  vhere 
one  of  the  factors  a^#  33/  or  33  could  have  a value  greater 
than  1.  It  is  also  evident  that  with  the  addition  of  a 
munber  of  secondary  arrivals  it  is  not  possible  to  deter- 
mine the  actual  value  of  a for  a particular  secondary 
arrival  from  a single  record. 
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Data  Sources 

Seismograms  of  selected  Nevada  Test  Site  events  were 
obtained  from  Lawrence  Livermore  Laboratories  and  Sandia 
Laboratories.  Seismograms  of  the  Massachusetts  Mountain 
earthquake^  vdnich  occurred  within  the  Nevada  Test  Site, 
were  obtained  from  Lawrence  Livermore  Laboratories 
(Rohrer  and  Springer,  1972;  Fis'ier  and  others,  1972). 
Records  of  other  earthquakes  occurring  within  or  near  the 
Ne'fada  Test  Site  and  thus  having  propagation  paths  similar 
tc  the  nuclear  events,  were  requested  but  were  not  avail- 
able. 

The  location  of  the  Nevada  Test  Site,  the  Lawrence 
Livermore  stations,  and  the  Sandia  stations  are  shown  in 
Figure  11.  ELK,  KNB,  LAC,  and  MNV  are  Lawrence  Livermore 
stations,  and  BMN,  DAC,  ELY,  LEE,  KEL,  and  TPH  are  Sandia 
stations.  Table  2 is  a listing  of  the  stations  and  their 
locations.  The  station  coverage  with  respect  to  Nevada 
Test  Site  is  illustrated  in  Figure  12. 

The  nuclear  explosion  records  obtained  for  analysis 
include  a multiple  event,  BLENTON/THISTLE , a single  event, 
DIDO  QUEEN,  and  five  events  A through  E of  unknown  nature. 
Information  relating  to  BLENTON/THISTLE  and  DIDO  QUEEN  is 
given  in  Table  3.  Table  4 is  a listing  of  the  depth  of 
burial  and  shot  medium  for  the  events  of  unknown  nature. 
Typical  records  of  the  events  analyzed  are  displayed  in 
Appendix  A. 
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Figure  11 


Location  of  Nevada  Teat  Site  and  Lawrence 
Livermore  (ELK,  KNB,  MNV,  and  LAC  and  Sandia 
(BMN,  DAC,  ELY,  LEE,  NEL,  and  TPH)  stations. 


DISTANCE  (KM) 
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Figure  12.  Station  coverage  of  Nevada  Test  Site. 

Distance  and  azimuth  stations  with  respect 
to  the  Nevada  Test  Site. 
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TABLE  2 

Seismic  Stations 


Code 

Station 

Latitude 

Longitude 

BMN 

Battle  Mountain/  Nev. 

40®25'53.3"N 

117‘'13‘18.4 

"W 

DAC  • 

Darwin#  California 

36®16'37.2"N 

117®35' 37.2 

"W 

ELK 

Elko#  Nevada 

40®44' 41.4"N 

115®14'19.6 

"W 

ELY 

Ely#  Nevada 

39®07'52.8"N 

114®53' 31.2 

"W 

KNB 

Kanab , Utah 

37®00'59.8"N 

112“49' 20.7 

"W 

LAC 

Landers#  California. 

34®23'23.2"N 

116®24'41.4 

"W 

LEE 

Leeds#  Utah 

37®14'34.8”N 

113®22' 36.0 

"V7 

MNV 

Mina#  Nevada 

38®25'56.0"N 

118®09'15.8 

"W 

NEL 

Nelson#  Nevada 

35®42'44.,0"N 

114®50'37.0 

"W 

TPH 

Tonopah#  Nevada 

■ 38®04'29.0”N 

117®13'21.0 

“W 
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TABLE  3 

BLENTON/THISTLE  and  DIDO  QUEEN  Data 


t i 


Date 

Shot  Tima 
Location 


Depth  of 
Burial 

Shot  Medium 


BLENTON 

4/30/69 

1700:00.04 

37“04' 53.4“N  . 
116®00'50.2“W 

• 56  km 


Tuff 


THISTLE 

4/30/69 

1700:00.04 

37®05' 25.0"N 
116®00’20.3"W 

.56  km 


Tuff 


DIDO  QUEEN 

6/5/73 

1700:00.00 

37“11' 06.1"N 
116®12'54.  l‘'W 

.39  km 


Tuff 


i;,'  !j 
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TABLE  4 


Depth  ar,d  Shot  Medium  for  Events  A-E 


Event 

Depth 

Shot  Medium 

■ A 

.24  km 

Alluvium 

B 

.56  km 

Tuff 

C 

.31  km 

Tuff 

D. 

.44  kmi 

Tuff 

£ 

.39  km 

Alluvium 
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Geologic  Setting 

Nevada,  southeastern  California,  and  western  Utah  lie 
within  the  Basin  and  Range  provinvoe  described  by  Thornbury 
(1965) . This  province  is  characterized  by  north-south 
trending  mountain  ranges  and  intermontane  basins  formed  as 
a result  of  block  faulting.  The  Basin  and  Range  province 
is  bounded  on  the  west  by  the  Sierra  Nevada  and  on  the 
east  by  the  Colorado  Plateau.  The  physiographic  provinces 
of  the  Nevada,  California,  and  Utah  area  are  outlined  in 
Figure  13. 

The  Nevada  Test  Site  is  located  in  a structurally 
complex,  highly  faulted  segment  of  the  Basin  and  Range 
province.  The  geologic  and  physiographic  features  of  the 
Nevada  Test  Site  and  surrounding  area  are  shown  in  Figurti 
14.  The  test  site  is  approximately  fifty  miles  east  of 
the  Cordilleran  eugeosyncline  eastern  boundary  and  about 
one  hundred  miles  west  of  the  miogeosyncline  eastern 
boundary.  The  area  is  also  north  and  east  of  the  Las 
Vegas  Valley-Walker  Lane  shear  zone. 

The  events  analyzed  in  this  study  were  detonated  in 
the  Yucca  Flats  area  of  the  Nevada  Test  Site,  and  the  epi- 
center of  the  Massachusetts  Mountain  earthquake  was  located 
in  the  southern  part  of  this  area.  A generalized  geologic 
map  of  the  Yucca  Flats  area  with  the  location  of  the 
earthquake  epicenter  and  the  area  in  which  the  nuclear 
events  were  detonated  are  given  in  Figure  15,  Upper 


•f  I 


Figure  14.  Index  map  of  the  geologic  and  phy 
graphic  features  of  the  Nevada  Te 
and  surrounding  area  (modified  fr 
Ekren,  1968,  p.  12) . 
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Figure  15.  Generalized  geologic  map  of  Yucca  Flats 
area  with  epicenter  of  Massachusetts  Mt. 
earthquake  and  area  in  which  nuclear  events 
were  detonated  (modified  from  Hinrichs, 
1968,  p.  242)  . 
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Procambrian  and  Paleozoic  sedimentary  rocks,  Mesozoic 
in'trusives.  Tertiary  volcanics,  and  Quaternary  alluvium 

t 

1 

are  exposed  in  the  area.  Quartzite  is  the  major  lithology 
I of  the  Upper  Precambrian  and  lower  part  of  the  Cambrian. 

I The  remaining  portion  of  .the  Paleozoic  is  composed  of 

I 

I carbonates,  conglomerates,  quartzites,  and  shales. 

Granitic  intrusives  are  the  only  Mesozoic  rocks  in  the 
area'.  Tertiary  rocks  consist  mainly  of  rhyolitic  tuffs. 
The  total  thickness  of  the  Yucca  Flats  stratigraphic  sec- 
tion is  44,500  .feet  and  is  divided  as  follows;  Pre- 
^ Cambrian  and  Paleozoic  rocks  38,800  feet;  Tertiary  rocks 

* 3,700  feet;  and  Quaternary  sediments  2,000  feet  (Hinrichs, 

1958) . 
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DIGITAL  ANALYSIS 

The  seismograms  were  prepared  for  computer  processing 
by  digitizing  the  records  using  a Computer  Equipment 
Corporation  Digitizer.  A FORTRAN  program  v/ritten  by 
'Taylor  (1974)  was  used  to  obtain  the  cepstra  and  a 
listing  is  provided  in  Appendix  B.  The  program  digitally 
performs  the  following  operations  on  the  input  data: 

1.  Interpolation  of  data  at  equal  increments. 

2.  Normalized  autocorrelation. 

3.  Calculation  of  Fourier  transform. 

4.  Calculation  of  amplitude  spectra. 

5.  Determination  of  log  spectra. 

6.  Normalized  autocorrelation. 

7.  Calculation  of  Fourier  transform. 

8.  Calculation  of  amplitude  spectra,  which  is  the 
cepstrum. 

An  example  of  the  results  of  the  different  operations 
performed  to  obtain  a cepstrum  is  shown  in  Figures  16 
through  19.  The  normalized  autocorrelation,  the  amplitude 
spectrum,  the  log  spectrum,  the  normalized  autocorrelation 
of  the  log  spectrum,  and  the  cepstrum  are  given  for  the 
BLENTON/THISTLE  event  recorded  at  the  KNB  station. 

The  input  data  were  interpolated  at  an  increment  of 
.05  second.  This  increment,  or  At,  was  adequate  to  satisfy 
the  Scimpling  Theorem 

At  < ^ (Huang,  1966)  (23) 
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where  is  the  maximum  frequency  present  in  the  signal. 
The  Nyquist  or  folding  frequency  is  then 


X 


N 


1 

2At 


10  Hz 


(24) 


or 

= 2TTfjj  = 62.8  radians/sec. 


(25) 


The  values  of  a Fourier  transform  computed  digitally 
differ  from  the  true  transform  values  by  a factor  of  At. 
With  the  application  of  two  transforms  to  the  data,  the 
computer  cepstrum  values  differ  from  the  true  cepstrum 
values  by  a factor  of  Aw,  the  frequency  increment.  The 
values  of  the  Fourier  transform  and  the  transform  of  the 
transform  were  digitally  computed  using  2048  points.  The 
value  of  the  frequency  increment  is 


Au 


2~b"4 '8 At  ~ •0614  radians/sec. 


(26) 


Denoting  the  true  cepstrum  value  as  C and  the  computer 
value  as  Cc, 


C = AuCc. 


(27) 


The  normalized  autocorrelation  which  was  performed  on  the 
data  introduces  a normalization  factor  into  the  computation 
for  the  true  cepstrum  value,  and  it  is  also  necessary  to 
take  the  square  root  of  the  computer  value  because  of  the 
property  that  the  Fourier  transform  of  an  autocorrelation 


—•Nj..  ■ • 


equals  the  square  of  the  amplitude.  A final  factor  to 
be  accounted  for  is  that  in  the  determination  of  the  log 
spectra,  the  logarithm  to  the  base  10  v/as  used,  but  the 
logarithm  expansion  of  equation  (9)  holds  for  the  logarithm 
to  the  base  e.  Thus,  the  computer  values  must  be  multi- 
plied by  a factor  of  2.3.  The  expression  for  the  true 
cepstrum  values  then  becomes 

C = 2.3Aw>^  (28) 
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The  factors  necessary  in  the  calculations  for  the  a 
scaling  factor  were  not  completely  understood  until  the 
majority  of  the  data  had  been  processed.  It  was  not 
possible  to  reprocess  data  to  obtain  the  normalization 
factors.  The  a factors  were  only  computed  for  the  arti- 
ficial double  events  created  using  the  DIDO  QUEEN  records. 


f 
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Results  and  Discussion 

It  has  been  observationally  established  that  signifi- 
cant energy  will  follow  the  arrival  of  any  P wave,  and 
this  following  energy  has  been  termed  the  P wave  coda. 
Secondary  arrivals  will  appear  within  the  coda  and  suffer 
from  the  coda's  contaminating  influence.  In  addition, 
the  digital  spectra  of  the  initial  P arrival  will  exhibit 
effects  due  to  contamination  of  the  coda.  Therefore,  in 
the  absence  of  any  well  defined  secondary  arrivals,  varia- 
tions in  the  cepstra  are  to  be  expected.  The  ambiguity 
in  interpretation  resulting  from  these  peaks  will  depend 
upon  the  degree  to  which  the  secondary  arrival  peaks  exceed 
the  extraneous  peaks.  Greenfield  (1969)  measured  coda 
power  levels  on  the  order  of  0. 1-1.0  of  the  primary  P 
levels.  The  figure  provided  by  Greenfield  actually  indi- 
cates power  level  ratios  on  the  order  of  0.01-0.1,  but  it 
appears  that  an  error  exists  in  the  labeling  of  the  figure. 
Ratios  on  the  order  of  0.1-1. 0 would  be  exactly  the  same 
order  as  that  expected  for  actual  secondary  arrivals. 

Thus,  the  identification  of  peaks  relating  to  secondary 
eirrivals  may  be  difficult,  and  may  prove  to  be  somewhat 
subjective. 

In  addition  to  the  problem  of  extraneous  peaks  it  is 
necessary  to  consider  the  problem  of  what  time  domain 
window  length  to  select  for  the  processing  of  the  data. 

A window  length  too  short  would  eliminate  desired  second 
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arrivals  / and  a window  length  too  long  v/ould  pass  addi- 
tional contaminating  energy  with  ro  increase  in  desired 
signal.  It  does  not  appear  that  any  theoretical  basis 
for  selection  of  an  optimum  window  length  exists.  Lacking 
a theoretical  basis#  the  effects  of  window  length  are 
investigated  in  the  following  section. 

Effects  of  Time  Domain  Window  Length 
The  combination  event  BLENTON/THISTLE  provides  an 
ideal  example  for  the  study  of  window  length  effects. 

As  previously  noted#  the  event  is  known  to  be  a simple 

« 

horizontal  dipole  and  the  records  of  the  event  should 
contain  a second  arrival  from  the  more  distant  event. 
Analysis  of  this  event  will  also  provide  an  indication 
of  the  magnitude  associated  with  cepstral  peaks. 

Cepstra  were  generated  for  8#  10#  15#  and  20  second 
windows  for  the  recordings  at  DAC  and  KNB . The  results 
are  shown  in  Figures  20  and  21.  The  peaks  of  the  DAC 
cepstra  are  consistent  for  the  different  window  lengths. 

The  peaks  of  the  KNB  cepstra  are  consistent  in  the  beginning 
two  seconds#  with  some  of  the  peaks  in  the  15  and  20 
second  windows  becoming  less  pronounced.  It  also  appears 
that  as  the  window  length  of  the  KNB  cepstra  is  increased 
more  peaks  occur  beyond  2 seconds. 

• On  the  basis  of  these  results  a window  length  of  J.0 
sec  was  selected  as  the  optimum  window.  A careful  study 
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Figure  20. 


Effect  of  window  length  on  cepstra  of 
BLENTON/THISTLE  event.  Window  lengths  of 
8,  10,  and  IS  seconds  for  signal  recorded 
at  KNB  and  DAC. 
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Figure  21.  Effect  of  window  length  on  cepstra  of 
BLENTON/THISTLE  event.  Window  lengths 
of  15  and  20  seconds  for  signal  recorded 
at  KNB  and  DAC. 
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of  Figure  20  v;ould  suggest  that  an  8 second  v/indow  is  more 
ne'arly  optimum.  It  is  not  expected  that  such  minor 
changes  in  windov;  length  can  greatly  affect  the  final 
results.  The  data  in  this  study  was,  therefore,  processed 
using  a window  length  of  10  seconds.  The  information 
contained  in  the  first  10  seconds  of  a signal  was  con- 
sidered  sufficient  since  the  secondary  arrivals  of  interest 
are  expected  to  occur  in  the  first  3 seconds  of  the 
record. 

In  addition  to  the  effects  of  window  length.  Figures 
20  and  21  indicate  the  problems  that  were  previously 
anticipated  in  the  identification  of  cepstral  peaks. 

While  the  relative  yields  of  the  two  events  is  not  known, 
the  expected  arrival  time  of  the  second  can  be  calculated 
and  is  indicated  by  the  vertical  arrows.  It  is  evident 
that  in  all  cases  a,  peak  does  exist  but  in  no  case  is  it 
a pronounced  peak. 

Results  from  a, Single  Event  (DIDO  QUEEN) 

and  an  Artificially  Created  Double  Event 

Cepstrura  analysis  was  applied  to  two  records  of  the 
single  event  DIDO  QUEEN.  There  were  not  enough  records 
of  processing  quality  available  to  analyze  delay-distance 
and  delay-azimuth  relationships.  The  cepstra  of  DIDO 
QUEEN  calculated  from  recordings  at  DAC  and  LEE  are  shown 
in  Figure  22.  Peaks  occur  in  the  DAC  capstrvim  at  .70, 


.95,  and  1.10  seconds  and  in  the  LEE  cepstrum  at  .75  and 
I.'IS  seconds.  The  peaks  occurring  at  .70  and  .75  seconds 
could  be  associated  v/ith  the  arrival  of  pP,  but  more 
records  are  needed  to  make  a firm  determination. 

The  DIDO  QUEEN  records  were  used  to  artificially 
create  the  record  of  a double  event.  An  a value  and  a 
delay  time  value  were  chosen,  the  digitized  signal  v/as 
multiplied  by  a,  delayed  by  the  time  value,  and  then 
added  to  the  original  digitized  input  to  artificially 
create  the  signal  for  a double  event.  Cepstra  of  the 
artificially  created  double  event  were  determined  for 
various  a and  delay  values.  The  LEE  cepstra  are  shown  in 
Figure  23  and  the  DAC  cepstra  in  Figure  24.  The  peaks  for 
the  artificially  created  second  direct  P occur  at  the 
chosen  delay  time  or  .05  seconds  later. 

The  a values  for  the  artificial  second  direct  P 
were  computed  using  equation  (31)  and  compared  with  the 
a value  used  in  creating  the  artificial  signal.  Table  5 
is  a listing  of  the  actual  a values  and  the  corapured  a 
values. 

It  is  evident  from  Table  5 that  the  formula  developed 
for  the  calculation  of  a is  valid,  although  the  percent 
error  varies  from  37%  high  to  32%  low.  A major  contribu- 
tion to  the  error  percentage  probably  results  from  the 
various  data  smoothing  operations  included  in  the  cepstral 
program.  It  does  not  appear  that  the  effects  of  these 
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Figure  23.  Cepstra  for  artificially  created  double 
event  using  DIDO  QUEEN  record  at  DAC. 
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Figure  24.  Cepstra  for  artificially  created  double 

, event  using  DIDO  QUEEN  record  at  DAC. 
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TADLE  5 

' Comparison  of  Actual  a Values  and  Computed  a Values 
for  Artificially  Created  Double  Events 

I • 

. I 

Delay 


Station 

Time 

(sec) 

Actual 

a 

Computed 
• a 

Reciprocal 

Root 

% Error 

DAC 

.50 

3.00 

3.61 

.277 

20 

DAC 

.50  • 

.333 

.263 

3.80 

21 

DAC 

1.00 

3.00  • 

2.87 

.349 

4 

DAC 

1.00 

1.00 

• 1.22 

.817 

22 

DAC 

1.00 

.333  . 

.226 

4.42 

32 

LEE 

.50 

3.00 

4.12 

.242 

22 

LEE 

.50 

.333 

,296 

3.38 

9 

LEE 

.90  , 

3.00 

3.66 

.273 

37 

LEE 

.90 

.33;3 

.363 

2.75 

11 

LEE  . 

1.00 

1.00 

.701 

1.42 

30 
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oporations  may  be  accounted  for  on  the  basis  of  any 
reasonable  theoretical  treatment.  To  eliminate  the 
smoothing  operations  results  in  an  increased  noise  level 
more  detrimental  to  calculations  of  a than  the  effects  of 
smoothing.  The  average  a calculated  on  the  basis  of 
several  stations  should  yield  a reasonable  estimate  of  a. 
The  a value  obtained  on  the  basis  of  a single  station  can 
be  considered  only  an  estimate  of  order.  As  previously 
noted,  it  is  also  evident  from  Table  5 that  without 
additional  information  as  to  the  nature  of  a secondary 
arrival  it  cannot  be  determined  which  root  is  the  actual 
value  of  a.  The  recovered  a values  may  therefore  be  of 
limited  value. 

The  peaks  occurring  in  the  cepstra  of  the  artificially 
created  series  are  an  indication  of  the  sensitivity  of  a 
cepstral  peak  to  an  a value.  The  peaks  having  a values 
as  low  as  0,333  were  clearly  visible  on  the  cepstra. 

Thus,  it  appears  that  secondary  arrivals  having  a values 
of  at  least  .333  have  the  potential  of  being  detected  by 
the  cepstrum  method. 

■fj 

, Results  from  a Known  Double  Event  (BLENTON/THISTLE) 

The  BLENTON/THISTLE  event,  which  was  previously  used 
in  a limited  manner  to  study  the  effects  of  window  length, 
provides  an  excellent  starting  point  for  any  attempted 
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analysis  of  multiple  arrivals-  The  explosions  consisted 
of  two  shots  which  were  detonated  simultaneously  and 
separated  by  1.22  km.  The  trend  of  the  line  connecting 
the  epicenters  was  N37°E.  Additional  knovm  inforraation 
concerning  the  shots  is  listed  in  Table  3.  On  the  basis 
of  the  known  information  it  would  be  expected  that  at 
least  three  significant  secondary  arrivals  would  exist  with 
delays  of  P^2“^l'  there  is  the 

potential  for  other  later  arrivals  resulting  from  slapdown 
effects.  The  P2”^l  should  exhibit  the  previo'tsly 

discussed  azimuthal  and  distance  dependences.  The  pP2"^l 
-P2_  delays  should  individually  exhibit  no  azimuth 
dependence,  but  may  yield  sectors  of  ambiguous  identifica- 
tion because  of  the  1.22  km  offset. 

Cepstrum  analysis  was  applied  to  nine  BLENTON/THISTLE 
records.  The  cepstrs  are  shown  in  Figures  25,  26,  and  27. 
The  delay  of  the  peak  of  greatest  amplitude  of  each  cepstra 
and  the  delay  of  other  peaks  of  lower  amplitude  were 
plotted  versus  the  distance  to  the  recording  station. 

When  two  or  more  peaks  have  the  same  ai'nplitude  and  are 
the  peaks  of  greatest  amplitude  occurring  on  the  cepstrum, 
.these  peaks  are  all  plotted  as  peaks  of  greatest  amplitude. 
The  delay-distance  plot  for  - BLENTON/THISTLE  is  shown  in 
Figure  28.  There  are  peaks  occurring  in  the  .40-. 50  sec 
range  on  the  cepstra  of  six  stations,  and  five  of  -tliese 
peaks  are  the  peaks  of  greatest  amplitude  for  that  cepstrxm. 
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, Figure  26.  BLENTON/THISTLE  cpi>3tra  for  I4NV,  NEL, 

! ELY,  and  LAC. 
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Figure  27.  BLENTON/THISTLE  cepstra  for  DAC,  KNB,  and  BMN 


TPH  DAC  NEL  LEE  ELY  KNB  LAC 


DISTANCE  (KM) 


-67- 


Also,  six  stations  have  ccpstral  peaks  in  .75-1.05  see 


ran'jc,  and  4 stations  have  a later  arrival  in  the  1.55- 


1.70  sec  range.  No  peaks  arc  visible  on  the  HNV  cepstrura. 


and  the  first  peaks  on  the  NEL  and  KND  cepstra  do  not 


occur  until  .75  and  .65  seconds,  respectively. 


Figure  29  is  the  delay-azimuth  plot  of  the  BLENTON/ 
THISTLE  secondary  arrivals.  The  trend  of  the  line  connec- 
ting the  midpoint  of  the  epicenters  and  the  theoretical 


curve  of  P2“P]l  for  a source  separation  of  1.22  km  is  also 
shown  on  the  plot.  kNV  lies  on  an  azimuth  perpendicular 


to  the  trend  of  the  line  connecting  the  midpoint  of  the 


epicenters,  TPH  and  NEL  lie  within  15°  of  the  perpendicu- 


lar. The  absence  of  arrivals  at  MNV  and  at  NEL  within 


.75  sec  is  consistent  with  the  expected  delay-azimuth 


relationship  for  two  horizontally  separated  events 


Using  the  .40-. 50  sec  arrivals  and  considering  the  absence 


of  arrivals  at  MNV  and  NEL,  the  delay-azimuth  relationship. 


for  two  horizontally  separated  events  is  visible  on  the 


plot.  The  delay-azimuth  relationships,  together  with  the 
delay-distance  relationship,  indicate  that  the  arrival  at 


.40-. 50  sec  is  the  second  direct  P arrival  from  the 


multiple  source. 


From  an  analysis  of  the  arrivals  in  the  .75-1.05  sec 


range,  an  increasing  delay  with  distance  relationship  is 
indicated.  This  arrival  then  may  be  the  pP  phase  from 


one  of  the  events.  The  later  arrival  at  1.55-1.70  sec 


could  possible  be  a Pg  phase,  but  a firm  determination 
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Figure  29.  Delay-azimuth  plot  for  BLENTON/THISTLE . 
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could  bo  accounted  for  by  either  of  these  fault  systems. 
FoXir  records  of  the  Massachusetts  Mountain  earthquake  v/cre 
available  for  analysis.  The  cepstra  are  shown  in  Figure 
30.  The  KNB  and  LAC  cepstra  show  peaks  with  amplitudes 
of  3 or  less,  while  the  MNV  and  ELK  cepstra  have  peaks 
with  cunplitudes  ranging  up  to  8.  The  delay-distance 
plot  in  Figure  31  and  the  delay-azimuth  plot  in  Figure  32 
have  a scatter  of  points,  and  it  is  not  possible  to  dis- 
tinguish a pattern  for  a secondary  arrival  on  the  basis 
of  only  four  stations. 

Despite  the  limited  data  available,  it  may  be  noted 
that  the  cepstra  peaks  appear  primarily  in  a 0.5  sec  to 
1.0  sec  window.  If  the  published  depth  estimate  of  4.6  km 
is  accepted  pP  should  appear  beyond  the  1 sec  window. 

It  thus  appears  reasonable  to  assume  that  the  cepstral 
peaks  are  not  the  result  of  pP  but  rather  result  from 
multiple  generation  on  the  fault  plane  or  differential 
path  effects  during  propagation. 

Results  from  EVENT  A 

Cepstrum  analysis  was  applied  to  seven  EVENT  A 
records.  The  cepstra  are  shown  in  Figures  33  and  34.  The 
delay-distance  plot  is  given  in  Figure  35.  The  peaks  of 
greatest  eimplitude  occur  over  a range  of  .35-. 90  seconds, 
with  5 stations  having  peaks  of  greatest  amplitude  in  a 
.70-. 90  sec  range.  There  is  no  later  arrival  cons.i.stently 
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cannot  be  made  on  the  basis  of  only  four  cepstra  exhibiting 
peaks  in  this  range. 

The  results  from  BLENTON/TIIISTLE  indicate  that 
cepstruin  analysis  is  a potential  method  of  detecting  the 
secondary  arrivals  for  an  unknown  event/  and  that  it  v;as 
possible  to  use  the  delay-distance  and  delay-azimuth  rela- 
tionships to  identify  secondary  arrivals.  However,  the 
delay-distance  and  delay-azimuth  plots  may  prove  very 
difficult  to  interpret. 
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Results  from  the  Massachusetts  Mountain  Earthquake 
The  Massachusetts  Mountain  earthquake  occurred  on 
August -5/  1971.  The  epicenter  was  located  at  latitude 
36®54.94'N  and  longitude  115®59.42'W,  and  is  shown  on 
the  Yucca  Flats  map  of  Figure  15.  The  focal  depth  of 
the  earthquake  was  4.6  km, ’and  the  magnitude  of  the  event 
was  3.5.  The  epicenter  was  south  of  the  Yucca  Flats  area 
of  the  Nevada'  Test  Site;  thus  the  earthquake  had  a great 
circle  propagation  path  nearly  identical  to  that  of  the 
nuclear  shots  detonated  in  that  region.  Fischer  and 
others  (1972)  determined  the  nodal  fault  plane  solution 
to  be  a plane  striking  N22®W  and  having  a vertical  dip  or 
a plane  striking  N68°E  and  having  a 86®  dip  to  the  north- 
west. The  main  earthquake  occurred  in  an  area  having  both 
northwest-trending  right-lateral  faults  cind  northeast- 
trending left-lateral  faults,  and  the  focal  plane  solution 
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Figure  32.  Delay-azimuth  plot  for  Massachusetts  Mountain 
earthquake . 
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Figure  33.  EVENT  A cepstra  for  MNV,  NEL 


Figure  34.  EVENT  A cepstra  for  ELY,  LAC, 
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i 

occurring  which  would  indicate  a clapdown  arrival,  nor  is 
there  an  absence  of  arrivals  at  any  station  v;hich  would 
Indicate  a multiple  event.  A suggestion  of  a consistent 
I secondary  arrival  does  exist,  and  is  indicated  by  the  line 

■'  dra\m  through  five  of  the  delay  points  in  Figure  35.  The 

^ slope  of  this  line  yields  a moveout  of  0.15  sec  per  100 

km  and  may  result  from  pP. 

The  delay-azimuth  plot  shown  in  Figure  36  does  not 
aid  in  distinguishing  a relationship  for  a secondary 
arrival  because  of  the  lack  of  station  coverage  in  the 
90“  and  270®  azimuth  area.  The  relationship  for  a single 
I - event  or  for  horizontally  separated  events  could  be  drawn 

I on  the  delay-azimuth  plot  for  the  data  available.  It  is 

not  possible  then  to  predict  the  nature  of  EVENT  A. 


- 
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Results  from  EVENT  B 

Ten  EVENT  B records  were  available  for  analysis.  The 
cepstra  are  shown  in  Figures  37,  38,  and  39,  and  the  delay- 
distance  plot  in  Figure  40.  Eight  of  the  cepstra  have 
peaks  occurring  in  the  .40-. 55  second  range,  and  a second 
concentration  of  peaks  occurs  on  eight  of  the  cepstra  in 
the  .70-. 90  second  range.  The  TPH  and  NEL  cepstra  have 
no  peaks  occurring  until  .80  and. 75  sec,  respectively,  and 
considering  this  absence  of  a .ivals  and  the  171®  dif- 
ference in  their  azimuth  with  respect  to  the  Nevada  Test 
Site,  TPH  and  NEL  may  possibly  be  located  on  azimuths 
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Figure  36.  Delay-azimuth  plot  for  EVENT  A. 
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pcrpcndicular  to  the  trend  of  a line  connecting  the  epi- 
center of  a double  event.  Considering  the  peaks  in  the 
.40-. 55  sec  range  of  the  delay-azimuth  plot  of  Figure  41, 
the  pattern  for  two  events  separated  horizontally  can  be 
drawn  on  the  plot.  The  concentration-  of  peaks  in  the 
.70-. 90  sec  range  may  be  associated  with  a pP  arrival. 

A line  possibly  related  to  pP  is  shown  on  Figure  40.  If 
this  line  is  correct,  it  indicates  a pP  movement  of  0.05 
sec  per  100  km.  There  is  also  a suggestion,  indicated  by 
the  dashed  line  in  Figure  40,  of  rn  event  arriving  v/ith 
relatively  constant  delay.  This  could  be  a slapdown 
signal.  v 

On  the  basis  of  the  delay-distance  and  delay-azimuth 
relationships  it  appears  that  EVENT  B is  a horizontally 
separated  dipole.  The  separation  appears  to  be  approxi- 
mately 1 km  and  trends  in  a N60®E  direction.  A slapdown 
signal  was  possibly  generated  1.23  seconds  following 
detonation. 


Results  from  EVENT  C 

Cepstrum  analysis  v.’as  applied  to  nine  EVENT  B 
records.  The  cepstra  are  shown  in  Figures  42  and  43  and 
Figure  44  is  the  delay-distance  plot  for  EVENT  B.  Seven 


of  the  cepstra  have  their  first  peaks  occurring  in  the 
.55-. 70  sec  range.  A line  through  these  points  may  be 
related  to  pP  and  is  shown  on  Figure  44.  If  correct. 


O Less  predominant  peaks 

Possible  delay-distance  relationship 

Possible  trend  of  epicenter  line 


Figure  41.  Delay -azimuth  plot  for  EVENT  B. 
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this  lino  indicates  a pP  movcout  of  0.07  see  per  100  km. 
The  two  other  cepstra,  LEE  and  ELK,  have  first  peaks 
occurring  at  .4,'S  sec.  None  of  the  stations  exhibit  an 
absence  of  peaks  which  v;ould  indicate  a multiple  event. 
The  delay-azimuth  plot  is  given  in  Figure  45.  It  appears 
that  EVENT  C is  a single  event  with  a pP  arrival  in  the 
.55  to  .70  sec  range.  DAC,  NEL,  ELY,  and  LAC  have  an 
arrival  in  the  1.5-1. 7 range  which  could  possible  be  a 
slapdown  cirrival. 


Results  from  EVENT  D 

Cepstrum  analysis  was  applied  to  nine  EVENT  D 
records.  The  cepstra  are  shown  in  Figures  46,  47,  and 
48,  and  the  delay-distance  plot  is  given  in  Figure, 49. 

Six  of  the  cepstra  have  first  peaks  occurring  in  the 
.45-. 60  second  range.  A line  of  very  small  slope  can  be 
drawn  on  the  delay-distance  plot  through  these  points, 
and  this  suggests  a multiple  event.  LAC  has  a first  pea.'*, 
at  .30  sec,  and  ELK  and  LEE  have  no  peaks  until  ,75  and 
.85  seconds,  respectively.  From  the  delay-azimuth  plot 

^ • 

in  Figure  50  it  is  evident  that  it  is  possible  for  LEE 
to  lie  on  an  azimuth  perpendicular  to  the  trend  of  a line 
coinccting  the  epicenters  of  a double  event,  and  the 
pattern  for  two  events  separated  horizontally  could  be 
drawn  through  the  data  point.  The  theoretical  curve  for 
horizontally  separated  events  having  a P-P  delay  of  .50 
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Figure  45.  Delay-azimuth  plot  for  EVENT  C. 
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Figure  47.  EVENT  D cepstra  for  ELK,  LEE,  and  KNB. 
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Figure  50.  Delay-azimuth  plot  for  EVENT  D. 
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:'.co  in  alno  :;hown  on  Fi'juro  50.  It  in  alno  por;:;ihlc  that 
there  is  a pP  arrival  in  the  .6 5-. 75  sec  ranye  v/ith  a 
ir.ovcout  of  .05  see  per  100  km.  The  corresponding  circle 
with  center  at  the  epicenter  could  be  drawn  on  the  delay- 
azimuth  plot.  Without  data  from  the  ?60*-310"  azimuth 
interval  the  nature  of  EVENT  D cannot  be  predicted. 

Results  from  EVENT  E 

Seven  EVENT  E records  were  availcible  for  analysis. 
The  cepstra  are  shown  in  Figures  51  and  52,  the  delay- 
distance  plot  in  Figure  53,  and  the  delay-azimuth  plot 
in  Figure  54.  Four  of  the  cepstra  have  first  peaks 
occurring  in  the  .45-. 55  sec  rcuige.  This  could  be  an 
indication  of  an  arrival.  There  is  also  an  indication 
of  a possible  pP  arrival  in  the  .80-1.00  sec  range  with 
moveout  of  .09  sec  per  100  km.  Lines  of  delay-distance 
dependence  have  been  placed  on  Figure  53,  but  are  dashed 
because  of  limited  data.  The  possible  delay-azimuth  pP 
relationship  has  been  drawn  on  Figure  54.  The  lack  of 
data  in  the  330 *-80*  azimuth  range  makes  impossible  any 
prediction  of  the  nature  of  EVENT  E. 


TPH  DAC  NEL  LEE 


-0-4' 


0-t4-0 


— 1 

m 

• 

<N 

r 

O 

• 

CM 

1.5  - 
1.0  - 

1 

m 

• 

o 

(oas)  AViaa 

1 

-! 

DISTANCE  (KM) 


I 

.5 

Delay  (aac) 


A Predominant  peak! 

O Leas  predominant  peaks 
Possible  delay-distance  relationship 


Floure  54.  Dalay-asimuth  plot  for  EVENT  E. 
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Suqgestions  for  Further  Research 

The  use  of  cepstrum  analysis  in  connection  with  the 
establishment  of  delay  dependence  on  azimuth  and  distance, 
as  proposed  in  this  study,  should  be  extended  to  include 
additional  natural  earthquakes.  This  would  provide  a more 
valid  basis  for  determining  the  ability  of  such  procedures 
to  differentiate  earthquakes  from  nuclear  events. 

Any  future  application  of  cepstrum  analysis  should 
involve  the  recovery  of  the  scale  factor  "a".  This  would 
be  particularly  important  in  attempts  to  identify  pP  at 
first  zone  distances.  The  "a"  factor  might  also  prove  of 
some  value  in  eliminating  extraneous  peaks.  The  utility 
of  "a"  in  this  connection,  however,  will  be  severely  limited 
by  its  ambiguous  nature. 

The  use  of  the  pseudo- autocorrelation,  suggested  by 
Bogert,  Healy,  and  Tukey  (1963) , should  be  investigated. 

The  pseudo- auto correlation  would  be  particularly  valuable 
in  the  identification  of  the  pP  arrival.  This  method 
also  has  the  potential  of  improving  the  detection  of  all 
arrivals  as  it  should  reduce  the  overall  noise  level. 
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APPENDIX  A 

TYPICAL  ANALOG  RECORDS 
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APPENDIX  B 


CXPSTRUM  PROGRAM  LISTING 


TfC  START  OF  DATA  WINDOW  “FASU^rD  FROM  STA^T  DF  DATA 
THI  EMD  or  DATA  WINDOW 
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